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^ ' Abstract. We characterize the error distribution of BATSE GRB locations by mod- 

U , ehng the distribution of separations between BATSE locations and IPN annuli. We 

' determine error model parameters by maximizing likelihood and rank the models by 

their Bayesian odds ratios. The best models have several systematic error terms. The 
simplest good model has a 1.9 degree systematic error with probability 73% and 5.4 de- 
grees with probability 27%. A more complex model adds a dependence on the datatype 
used to derive the location. 
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: ANALYSIS METHOD 

p • The purpose of this paper is to develop better models of the error distribution 

^ \ of burst locations in the 3B [5] and 4B [6] catalogs, all of which were obtained 

^ I with the same algorithm. A more accurate error distribution will aid counterpart 

^ ■ searches and improve the accuracy of limits on burst repetition and clustering. 

Measuring the actual location errors requires a comparison dataset of sufficient 
accuracy and size — the most suitable set is the locations of the Interplanetary 
^ \ Network (IPN) [4]. The advantages of this dataset are its large size and small 

■ errors. The disadvantages are that most of the error boxes are single annuli and 

that the sample is biased towards bright GRBs. 

After removing 11 wide annuli (Scr error > 0.8°), the IPN supplement [4] to the 
4B catalog consists of 442 GRBs. Of these, 30 have intersecting annuli that yield 
essentially 'point' locations and thus provide measurements of the angles 7 between 
the true locations and the BATSE locations. For the remaining events with only 
single annuli we can determine only the closest approach angles p between the 
annuli and the corresponding BATSE locations. 

Our method of analyzing this set of 7 and p measurements is similar to that 
of Graziani & Lamb [3]. We assume various models for the systematic error a^js 
and obtain the total error atot as the rms sum of the statistical error a"stat listed in 
the catalogs and a^ys- For the probability distribution P(7) we assume the Fisher 
distribution, which is generally considered to be the equivalent of the Gaussian 



distribution on the sphere [2]. For those GRBs with single annuh for which we 
know p but not 7, we derive the exact distribution P{p) from P(7) [1]. 

The model parameters arc obtained by maximizing the likelihood, which is the 
product of the probability, according to the model, of the observations: 

* i 

Once the likelihood is optimized for each model, the models are compared by 
their Bayesian odds ratio [8]: 



O = = ^prior(B) X L(B) X Occam(B) 

P(A) Pprior(A) X L(A) X Occam(A) ^ ' 



We set the first factor, representing our prior preferences for the models, to unity. 
The second factor, the likelihood ratio, indicates how strongly the data favor one 
model over another. The "Occam's factor" penalizes a model to the extent that the 
better fit is obtained by additional parameters. For the data and models presented 
below, the factor per parameter by which a model is penalized ranges from about 
5 to 20 so that the dominant factors in distinguishing between the models are the 
likelihood ratios. 



RESULTS 

We first try the minimal model (Model 0), which has (Tgys = 1.6° for all GRBs. 

The value 1.6° was obtained using 50 BATSE locations and the corresponding 
WATCH, COMPTEL and IPN locations as the rms difference between the actual 
separations 7 and the errors estimated from the errors of the other instrument 
and the BATSE statistical errors [5,7]. Figure 1 shows that the agreement of the 
minimal model with the data is very poor. Model 1 generalizes Model by making 
(Tgys a free parameter with the same value for all GRBs. While a vast improvement 
over Model (sec Table 1), the agreement between Model 1 and the data is still 
quite poor (histogram not shown). 

While the minimal model is the previously published quantitative model, the 
poor agreement was expected by the BATSE Team, e.g., the 3B catalog [5] notes 
"However, the various tests yield only an estimate of the average systematic error. 
A small fraction of the locations could be substantially worse than the average, 
i.e., the location error distribution may have a non-Gaussian tail." We now have 
enough accurate comparison locations to test models with additional parameters. 

Model 2 uses the sum of two Fisher distributions with differing values of agys to 
implement a Gaussian distribution with an extended tail. It is a very major im- 
provement over the simpler models (Table 1) . Model 2 is an uncorrelated model — it 
does not use any datum to determine the value of agys for a location and is best 
thought of as a Gaussian with wings rather than a two-component model because 
BATSE data does not assign a burst to a component. 
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FIGURE 1. Histograms of the' data and models. Left: 30 events for which two anmih are 
available. Right: 412 events for which only single annuli are available. The data and both models 
are binned in units of fJtot.O) as calculated with the minimal model (Model 0) of (jgys = 1-6°. 
The solid lines show the data, the dotted lines the minimal model, and the dashed lines show 
Model 12. The histogram for Model 2 (not shown) is very similar to the histogram of Model 12. 



TABLE 1. The Models 

Odds 
Ratios 



Model Description Logio Parameter 

Likelihood Values 



Minimal 



1 

4 



One cr. 



sys 

^(astat/l°)" 



2 Core plus Tail 



12 Datatype 

Dependence 



140.8 
215.0 
217.4 

226.7 

229.9 



(Tsys = 1.6° fixed 

o-sys = 3.07° ± 0.12° 

A = 3.22° ±0.13° 
a = 0.18 ±0.05 

/i = 0.73 ± 0.08 

a,ys,i = l-9°±0.2° 
CTsys,2 = 5.4° ± 0.7° 



/coNT = 0.76 ±0.07 
4oNT = l-6°±0.2° 
^CONT = 5.4° ± 0.8° 

Mother = 3.7° ± 0.3° 
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1 X 10**3 
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6 X 10^° 
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2 X 10^2 
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1 X 10*5 



So far, the only significant correlation we have found between location errors 
and other location properties is with the datatype used to locate the event. Fit- 
ting Model 1, one value of CTgys, to locations obtained with the 16-channel CONT 
datatype, the best fit value of (Tgys is 2.81° ±0.14°, while for the other datatypes, the 
best value is 3.69° ±0.29°, clearly showing that locations obtained using 16-channel 
data are more accurate than those obtained with 4- or 1-channel data. 

Based upon this correlation. Model 12 was constructed: events located with 
CONT data have a core-plus-tail distribution similar to that of Model 2, while 
events located using other datatypes have a single value for Csys (Table 1). The 
data and Model 12 agree well (Fig. 1). The evidence for an extended tail for the 
"other" locations is marginal (odds ratio of 4). So far extensions of Model 12, such 
as adding intensity dependences, have at best marginal odds ratio improvements 
(< 10) and frequently unconvincing parameter values. 

We find only weak evidence for a direct dependence of (jgys on the intensity of a 
burst, as measured either by crgtat or by fluence. Model 4, with agys a power law of 
o'stat [3], is a definite improvement over Model 1. However, an uncorrelated model 
is yet superior, indicating that the intensity correlation of Model 4 is not supported 
by the current data. Additionally, a histogram (not shown) shows Model 4 to be 
a poor fit to the data. Consequently, the extrapolation of Model 4 to faint bursts 
(e.g., (Tsys = 8° for (Tstat = 10°) [3], and the consequent conclusion that the 3B 
catalog has larger systematic errors for faint bursts than the IB catalog [3] are not 
justified. 

While we do not find an intrinsic correlation of Usys with burst intensity, the 
fraction of events located with CONT data decreases with decreasing intensity, 
thereby causing an indirect intensity correlation via the datatype correlation. The 
fraction of CONT locations falls from 70% for events with CTgtat < 3.7° to 40% for 
the remaining events. However, this indirect correlation does not cause cTgys to scale 
with intensity, instead it causes the fraction of events with agys = 3.7° to be higher 
for sets of fainter bursts. 

While the odds ratio favors Model 12 over Model 2, Model 2 is probably sufficient 
for most purposes. Here we give a detailed specification of Model 2 so that others 
can implement these two models. 

Model 2 consists of two terms: 

P = /iPi ± (1 - /i)P2 (3) 

Each term is a Fisher probabihty distribution so that the probabihty of the true 
location lying in a ring of 71 to 72 about the BATSE location is: 

Pi = f [ ^^^e«^-«^ = rrf7sin7e'^^-^^ (4) 

By definition, the radius o"tot,i contains 68% of the probability. Setting Pi — 0.68, 
7i = 0, and 72 = Utot,?, we obtain (with Utot,? in radians): 



For each burst, the values of the total location errors crtot,i are calculated from 
the statistical error Ustat listed in the catalog and the model errors asys,*: 



CONCLUSIONS 

The best models are 2 and 12, which are implementations of the BATSE Team's 
expectation that the location error distribution is Gaussian modified with an ex- 
tended tail. Model 12 incorporates the only significant correlation that we have 
identified, basing the error distribution on the datatype used to derive the loca- 
tion. 

While the fraction of locations belonging to the larger systematic error term 
is 27% according to Model 2, this does not mean that 27% of the locations have 
large errors — the large systematic error term has a significant probability of a small 
location error. For events with (Jstat ^ crsys,i, the fraction of locations with an error 
past 2(Tsys,i = 3.8° is 16%. 

The results herein are upper-limits on BATSE systematic errors since they are 
based on the assumption that there are no errors in the comparison dataset. As part 
of this effort, most of the outliers were re-examined and eight significant revisions 
were made to the preliminary IPN catalog. Even a single outlier past 4(7 can 
significantly effect the error model parameter values. 

We continue to test additional models. Post-4B improvements in the location 
algorithm will lead to revised error model parameters. 
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